Renormalization approach to many-particle systems 



(N 
O 

o 

> 
O 



K. W. Becker, A. Hiibsch, and T. Sommcr 
Institut fur Theoretische Physik, Technische Universitdt Dresden, D-01062 Dresden, Germany 

(February 1, 2008) 

This paper presents a renormalization approach to many-particle systems. By starting from a 
bare Hamiltonian TL = Ho + Hi with an unperturbed part Ho and a perturbation Hi , we define 
an effective Hamiltonian which has a band-diagonal shape with respect to the eigenbasis of Ho- 
This means that all transition matrix elements are suppressed which have energy differences larger 
than a given cutoff A that is smaller than the cutoff A of the original Hamiltonian. This property 
resembles a recent flow equation approach on the basis of continuous unitary transformations. For 
demonstration of the method we discuss an exact solvable model, as well as the Anderson-lattice 
model where the well-known quasiparticle behavior of heavy fermions is derived. 
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I. INTRODUCTION 



If one considers a many-particle problem one would like 
to diagonalize the Hamiltonian so that observables of in- 
terest can be calculated. However, only a few models can 
be solved explicitly so that often one has to use approxi- 
mations or numerical approaches. Recently Wegner and 
coworkers [1,2] and Glazek and Wilson [3,4] proposed to 
solve the eigenvalue problem of a many-body system by 
transforming the Hamiltonian in such a way that it be- 
comes more and more diagonal. The method is based on 
the introduction of continuous unitary transformations 
and is formulated in terms of differential equations for 
the parameters of the Hamiltonian. These equations are 
called flow equations. 

In this paper we present a renormalization approach 
which is directly based on perturbation theory for the 
Hamiltonian instead of using a formulation in terms of 
differential equations. The starting point is a many- 
particle Hamiltonian H = Ho + Hi- It is assumed that 
the eigenstates and eigenenergies of the unperturbed part 
Ho are known. The interaction Hi usually prevents the 
solution of the full eigenvalue problem of H. In Sec. II, 
we shall derive an effective Hamiltonian H\ on the ba- 
sis of a perturbational approach. The obtained Hamil- 
tonian H\ has no matrix elements (with respect to the 
unperturbed Hamiltonian Ho) which belong to transi- 
tions larger than a given energy cutoff A. In contrast 
to standard perturbation theory no vanishing energy de- 
nominators appear. In Sec. Ill, the result for H\ will be 
used to formulate a A-dependent renormalization scheme 
by repeatedly applying a stepwise transformation on the 
effective Hamiltonian. In an infinitesimal formulation the 
renormalization of the Hamiltonian will be described by 
flow equations. At the end of this section a comparison 
of our approach with those of Wegner [1] and of Glatzck 
and Wilson [3,4] is given. For demonstration of the renor- 
malization approach in Sec. IV the exactly solvable Fano- 
Andcrson model [6,7] is discussed. This model describes 



two types of interacting electrons without correlations 
and the exact result will be obtained in the framework 
of the renormalization approach. Finally, in Sec. V the 
application to the Anderson-lattice model [6] is discussed 
in order to derive the well-established quasiparticle be- 
havior of heavy fermions. The conclusions in Sec. VI 
conclude the paper. 



II. EFFECTIVE HAMILTONIAN AND 
PERTURBATION THEORY 

Let us start from a decomposition of a given many- 
particle Hamiltonian H into an unperturbed part Ho and 
into a perturbation Hi 



H = H + eHi =:H{e). 



(1) 



We assume that the eigenvalue problem of Ho is known 

H \n)=E^\n). (2) 

H\ is the interaction. Its presence usually prevents a so- 
lution of the eigenvalue problem of the full Hamiltonian. 
The parameter e accounts for the order of perturbation 
processes discussed below. Let us define a projection op- 
erator P\ by 



PxA = 



E 

?(0) pCO) 



|n)(m| (n|A|m) 



(3) 



E^'\<\ 



Here, Pa is a superoperator acting on usual operators 
A of the unitary space. It projects on that part of A 
which is formed by all dyads |n)(m| with energy differ- 
ences l-E^ ' — Em | less or equal to a given cutoff A, where 
A is smaller than the cutoff A of the original model. Note 
that in Eq. (3) neither |n) nor \m) have to be low-energy 
eigenstates of Hq. However, their energy difference has 
to be restricted to values < A. Furthermore, it is useful 
to define the projection operator 



Qa 



(4) 



1 



which is orthogonal to P x . Q x projects on high-energy 
transitions larger than the cutoff A. 

The goal of the present method is to transform the ini- 
tial Hamiltonian Tt (with a large energy cutoff A) into an 
effective Hamiltonian TL X which has no matrix elements 
belonging to transitions larger than A. This is achieved 
by a unitary transformation so that the effective Hamil- 
tonian will have the same eigenspectrum as the original 
Hamiltonian TL. However, as it will turn out, it is espe- 
cially suitable to describe the low-energy excitations of 
the system. TL X is defined by 



He 



(5) 



and similar expression in higher order in e. Here Lo de- 
notes the unperturbed Liouville operator which is defined 
by lioA = [Hq, A] for any operator variable A. Note that 
by Eqs. (10), (11), and the equivalent equations for the 
higher orders in e the generators are not completely 

fixed. The reason is that the parts P^Xj-"' of X^ 1 ' are 
not fixed by Eq. (6) and can still be chosen arbitrarily. 
In the following we use for convenience 



P X X 



(l) 



^ A 



P X X 



(2) 



^ A 



0. 



(12) 



Inserting Eqs. (10) and (11) into Eq. (9) one finds for 
the effective Hamiltonian TL X up to second order in TL\ 



The generator X\ of the transformation has to be anti- 
Hermitian, X^ = —X x , so that TL X is Hermitian for any 
A. We look for an appropriate generator X\ so that H\ 
has no matrix elements belonging to transitions larger 
than A. This means that the following condition 



QaHa = 



(6) 



has to be fulfilled. Equation (6) will be used below to 
specify X x . 

The expression (5) for the effective Hamiltonian 7i x 
can we expanded with respect to X\ 



Hx = H + {X x ,H] 
1 



[Xx, [X\, H] 



3! 



[X x ,[X x ,[X x ,n]]] 



(7) 



By assuming that X x can be written as a power series in 
the perturbation parameter e 



X X = eX^ 



L X 



(2) 



S 3 X[ 3) 



(8) 



the effective Hamiltonian TL X can be expanded in a power 
series in e as well 



Hx=n + e{Hi+[x£ ) ,Ho\} 



(9) 



{ 



1 

+ 2! 



X 



(i) 

A ' 



X x 2 \n 
X x , Ho 



0(e 3 ). 



H x =Ho + PxH 1 -^P x 



(QaWi),=?-(QaWi) 

1^0 



(13) 



(PaWi),=^(QaHi) 

1^0 



where e was set equal to 1. Note that the perturbation 
expansion (13) can easily be extended to higher orders 
in e. One should also note that Eq. (13) automatically 
guarantees the correct size dependence of the Hamilto- 
nian due to the commutators appearing in Eq. (13). 

In the following section we will use Eq. (13) to estab- 
lish a renormalization approach to many-particle systems 
by successively reducing the cutoff energy A to smaller 
and smaller values. In the limit A — > expression (13) 
provides an expression for an effective Hamiltonian that 
acts on a possibly degenerate ground-state of the unper- 
turbed Hamiltonian. For this case one can also show that 
Eq. (13) reduces to an expression derived from usual per- 
turbation theory [8]. A systematic extension of Eq. (13) 
to higher-order perturbation theory can also be used to 
provide a simple scheme for an algebraic evaluation of 
higher terms of TLi X= o) by use of a computer [9]. Such a 
scheme does not require special properties of the eigen- 
value spectrum of the unperturbed Hamiltonian TLo- In 
particular, an equidistant spectrum is not needed. This 
property was implied for a perturbation expansion [10] 
based on Wegner's flow equation method. 



III. RENORMALIZATION APPROACH 



(n) 

The different contributions X^ to the generator of the 
unitary transformation (5) can be successively deter- 
mined from Eq. (6). One finds 



QaAT = 



i 1} = r (QaHx), 

-Li 



2L 



-^Qa 



[QxHi),— (QaHi) 
(PaWi) ,^(QaW x ) 

^0 



(10) 



(11) 



In this section we will discuss the elimination proce- 
dure for the interaction 7i\. However, instead of trans- 
forming the Hamiltonian in one step as was done in the 
preceding section the transformation will be performed 
successively. Or more formally spoken, instead of ap- 
plying the elimination of high-energy excitations in one 
step a sequence of stepwise transformations is used in 
order to obtain an effectively diagonal model. In an in- 
finitesimal formulation the renormalization of the cou- 
pling constants will be described by flow equations. In 
order to find these equations some approximations will 



2 



be necessary. The goal of the method is to transform 
the initial Hamiltonian (with a large cutoff A) into an ef- 
fective Hamiltonian which possesses only coupling terms 
between states with an energy difference less or equal to 
A. Thereby, the method yields rcnormalization equations 
as function of the cutoff A. 

Let us start from the renormalizcd Hamiltonian 



H\ = Hq.x + Hi, 



(14) 



after all excitations with energy differences larger than 
A have been eliminated. In the next step an additional 
renormalization of H\ to a new Hamiltonian 7Y(a-aa) 
is done by eliminating all excitations within the energy 
shell between the cutoff A and a somewhat smaller energy 
cutoff (A — AA). The new Hamiltonian 7Y(a-aa) will be 
calculated by use of the perturbation theory discussed 
above. One finds 



H(\- 



(A — AA) 



(15) 



' (\ — A\)H\ 

1 

-7TMA-AA) 



(A-AA) 



(Q(A-AA)^1,a), 
(P(A-AA)Wi,a), 



Lo,a 

1 



J 0,A 



(Q(A-AA)^1,a) 
(Q(A-AA)Wi,a) 



L ,a is now the Liouville operator with respect to the A- 
dependent unperturbed Hamiltonian Hq,\. As is obvious 
from Eq. (15) the second-order term on the right-hand 
side (rhs) may be divided into two parts. The first one 
connects eigenstates of Ho,\ with the same energy. This 
part commutes with Ho,\ and can therefore be considered 
as renormalization of the unperturbed Hamiltonian 



W , (a-aa) — Ho,\ — 
~o P o (Q(a-aa)^i.a) 



(16) 



J 0,A 



(Q(A-AA)^1,a) 



In contrast the second part connects eigenstates of Ho.x 
with different energies. It therefore represents a renor- 
malization of the interaction part of the Hamiltonian 



Hi 



(A-AA) 



(A-AA) 



(17) 



(A-AA) 



(P(A-AA)^1,a), 



J 0,A 



(Q 



(A-AA) 



Hi^x) 



4 P 



(A-AA) 



(Q(A-AA)^1,a), -z (Q(A-AA)^1,a) 



-*0,A 



(A-AA) 



(P(A-AA)Wi,a), j (Q(a-aa)^i,a) 



J 0, A 



We assume that AA is small. Therefore, only the mixed 
term, i.e., the second part on the rhs of Eq. (17), con- 
tributes to the renormalization of Hi,\. Thus, Eqs. (16) 



and (17) describe the renormalization of the unperturbed 
and of the interaction part of the Hamiltonian. They rep- 
resent the main result of our theoretical formalism. In the 
limit A A — > 0, i.e. for vanishing shell width, these equa- 
tion can be used to derive differential equations for the 
dependence of the Hamiltonian on the cutoff energy A. 
The resulting equations for the parameters of the Hamil- 
tonian are called flow equations. Their solution depend 
on the initial values of the parameters of the Hamiltonian 
and determine to final Hamiltonian in the limit A — > 0. 
Note that for A — > the resulting Hamiltonian only con- 
sists of the unperturbed part Hqi\— «-o)i so that an effec- 
tively diagonal Hamiltonian is obtained. The interaction 
Wwa-^o) completely vanishes since it is completely used 
up in the renormalization procedure. The detailed struc- 
ture of the flow equations strongly depends on the system 
under consideration. However, one can easily write down 
these equations in matrix representation with respect to 
the unperturbed Hamiltonian (compare Appendix A). 

Next we discuss the relation between the present renor- 
malization approach and Wegner's flow equation method 
[1]. The approach presented here is also based on the 
idea to formulate flow equations for the Hamiltonian. 
However, there are substantial differences between both 
methods: i) The present approach starts directly from 
the renormalization of the Hamiltonian and not from 
continuous unitary transformations in differential form, 
ii) Moreover, the band diagonal form of the transformed 
Hamiltonian results from a physically motivated projec- 
tion operator on low-energy transitions instead of intro- 
ducing cutoff functions to restrict the Hamiltonian ma- 
trix, iii) The result of the renormalization approach (16) 
and (17) is directly formulated in terms of an operator 
relation. Wegner's flow equation approach [1] is usually 
written in matrix representation though it is also possi- 
ble to give a formulation in terms of operators [5]. iv) 
As was discussed in Sec. II, there is a direct connection 
between the renormalization method and usual pertur- 
bation theory. 



IV. THE EXACT SOLVABLE FANO-ANDERSON 
MODEL 

In this section we illustrate the renormalization ap- 
proach for the case of a simple model, i.e., we apply the 
method to the exact solvable Fano- Anderson model [6,7]. 
This model is given by 



H = Hq + Hi, 



(18) 



i.rn k,rn 
= { £ f /km/km + £k 4m c km) , 



k , rn 



3 



Wl = ^ E Vk(/L Ckm e ikR -+h, 

k,i,m 



k,?n 

A possible realization of the model is a periodic system of 
localized / electrons which interact with conducting elec- 
trons thereby neglecting correlation effects. The index i 
denotes the / sites, k is the wave vector, and Vk describes 
the hybridization between conduction and localized elec- 
trons. The excitation energies £k and £/ for conduction 
and localized electrons are measured from the chemical 
potential fi. As a further simplification, both types of 
electrons are assumed to have the same angular momen- 
tum index m with values m = 1 . . .Vf. Of course, the 
model is easily solved and leads to two hybridized bands 



H 



E(«) t 

k,m 



km 



k.rn 



(19) 



Jk ~ 2 



,(<*>/?) 



4|V k | ! 



with eigenmodes aj CT)i and /3 km given by certain linear 
combinations of c krn and /^ m , i.e., 



a k.m ~ Wk / km + 



T " kC km' 4m 



^k/l 



km 



it k c 



km ' 



-k,A 



and 



£ k a wm depend on the energy cutoff A and 
moreover e£ x depend on the wave vector k. Next, we 
evaluate the new Hamiltonian Ho^a-aa) by projecting 
out all excitations between A — AA and A [compare (16)]. 
By inserting Eq. (21) into Eq. (16) we obtain 




^O.(A-AA) — T~lo,\ 
\Vk\ 2 , 



(22) 



/km/ 1 



km 



,Ckn 



where the Fermion commutation relations were used. 
The prime / above the sum indicates that the condition 
(A - AA) < |e k A - elj < A has to be fulfilled. Note 
that in Eq. (22) only second-order contributions are ex- 
plicitly given. Higher-order terms are indicated by dots 
.... The renormalization equations for the one-particle 
energies are easily found from Eq. (22) and read 



-/ 

'k,(A — AA) 



\Vkf 



-k,X 



-k,\ £ k,X 
2 



£ k.(A-AA) _ £ k,X 



where 



\Vk? 



£■ 



k,A 



'k.A 



6(A,AA), 

(23) 

1 6(A,AA) 
J (24) 



u k 



Vk 



1 1 

2 ~ 2~' 

1 1 



( £k - £/ ) 2 + 4|V k p 
Ek-£f 



(20) 



In the renormalization approach we are integrating out 
particle-hole excitations of conduction and / electrons 
described by the hybridization Jiy. We expect to obtain 
finally an effectively free model. 

The starting point of the method is a renormalized 
Hamiltonian TL\ which is obtained after all excitations 
with energies larger than a given cutoff A have been elim- 
inated. Due to the result of the preceeding section it 
should have the following form 



7i\ = H, 



0,A 



H 



1,A, 



Ho.X = J2 ( £ k,A /km/km + 4,A 4m < *' 
k,m 

Hl,X = E ^ k (/km C km + C k m/krn) 



(21) 



= Pa Wi. 



k , m 

|^.A-<a|<^ 



As it turns out no renormalization of the hybridization 
Vk occurs so that Vk is assumed to be A independent 
from the beginning. However, the one-particle energies 



6(A, AA) = 9(| 



£ k.x - £ 



k.A 



- (A - AA) 



-9 



£ Lx ~ £ k,x\ ~ A 



(25) 



To solve these equations we subtract Eq. (23) from Eq. 
(24), multiply both sides by (e k A — e£ A ), and divide by 
A A. By performing the limit A A — > we obtain 



d 

r/A 



(4,a-4,a) : 



AA^O AA 



Note that Eq. (26) includes all renormalization contri- 
butions of higher orders in the interaction. It turns out 
that the solution of Eq. (26) shows a step like behavior 
and reads 



l e k,A 
with 



-kjj 



lek-e/l+Oceflek-e/l-A) (27) 



Ck = J(ek-e f r + 4|y k |2-|e k -e/|. 



To prove this result one has to insert Eq. (27) into Eq. 
(26) which gives — 4|Vk| 2 d(|e k - e f \ - A) for both sides. 
Note that the evaluation of the rhs. of Eq. (26) together 



with Eq. (25) does not give -4| V k \ 2 S{\e c k x - e 
as would be expected of first glance, 
the step like renormalization of |e k A — 
Eq. (27). 

For A — > one finds 



/ 

k,Al 



A) 



This follows from 
£ L x I > described by 



4 



(£ k -e/) 2 + 4|y k | 2 



(28) 



where e k and £j£ denote the one-particle energies at 
A -> 0, 



ew = e 



k,(A->0)> 



-/ _ / 

£ k — £ k,(A-+0)' 



(29) 



On the other hand, it is easily seen that the summation 
of Eqs. (23) and (24) leads to 



si 



Sk + Sf. 



(30) 



The system of Eqs. (28) and (30) can easily be solved 
and gives 



-f _ £k + £/ _ 1 



-^( £k - £f y+4\v k \< 



(31) 



(32) 



Note that Eq. (32) agrees with the exact solution Eq. 
(19). As a result we have obtained the following free 
model 



(33) 



Note that the interaction Hi, a vanishes for A — * since 
it is completely used to renormalize the one-particle en- 
ergies e k and e£. The result Eqs. (33) and (32) can 
be used to evaluate expectation values. Since 7i(.\->o) 
emerged from the original model 7i by an unitary trans- 
formation, the free energie can also be calculated from 

H(A-+0) 



F = -IlnTre"^ = -- lnTre _/3W '^°> 

(3 /3 



Ft 



(A-0)- 



(34) 

Then, the / electron occupation number is found from 



F 



(A-0) 



by functional derivative 



1 OF 



(ft f . ) _ ±^f_ _ 



1 OF, 



(A^O) 



N def N def 
1 /l 1 



(35) 



£k 



1 + e^t V 2 



1 + e^Z 



2v/( £ k-£/) 2 + 4|^ k |V 

1 _ 1 £k - £/ 

2 2 v /( £k - e/ )2+4|F k | s 



where Eqs. (32) and (34) was used. On the other hand, 
static and dynamic expectation values can be evaluated 
by applying the unitary transformations not only on the 
Hamiltonian but on all operator quantities which ap- 
pear in expectation values. By exploiting that operator 
expressions under a trace are invariant against unitary 



transformations the / electron occupation number is also 
given by 

(/L/km> = (/LW/k™(A)>A ( 36 ) 
= (./ , L(A-0)/ km (A^0)) (A _ o) 



where 



fLW=^ ft, 



(37) 



Here, (...) A means the expectation value formed with Tt\ . 
One possible way to proceed is to derive flow equations 
for the A-dependent operators / km (A), / km (A). However, 
this may lead to inconsistent approximations which are 
necessary to solve the obtained flow equations. On the 
other hand, this problem can be avoided by using an ap- 
propriate ansatz for the A-dcpcndence. Therefore, for the 
A-dependent operators / km (A) and for the second opera- 
tor c km (A) we shall use 

/L(A)=«k,A/L + «k,AcL and (38) 

C L( A ) = ~ V K\ /km + M k,A 4m 

with |w k ,A| 2 + |-uic,A | 2 = 1. This form is suggested by re- 
peatedly applying the unitary transformation on / km and 
c km . With Eq. (38) the expectation value (/ km /km) can 
easily be evaluated. One finds 

</L/km) = (/L(A - 0)/ km (A - 0)) M (39) 

= K, ( a^)I 2 — ^7 + kk.cA^o)! 2 ^^. 



If one compares both results for (/ km / km ) the prefactors 
u k,(A-»o) ancl u k,(A^o) can be found 



|«k,(A^o)| 2 = i + ^ 

K,(A^0)| 2 



£ k - Sf 



1 1 

2 + 2" v /( £k - £/ )2 + 4 |y k p : 

11 £k - £/ 



(40) 



2 2^( £k _ £/ )2 + 4 |y k |2 



This result agrees with Eq. (20). Thus we have obtained 
the exact transformations for the / and c operators. It 
follows that all static and and dynamic quantities involv- 
ing electron creation and annihilation operators can be 
calculated in the framework of the renormalization ap- 
proach. 
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APPLICATION TO PERIODIC ANDERSON 
MODEL 



In the following section we will refer to the Anderson 
lattice model [6] which describes an interacting system 
of conduction electrons and correlated localized 4/ elec- 
trons, arranged periodically on a lattice. Within a simpli- 
fied version the Hamiltonian of the model can be written 
as 

H = H +Hi, (41) 

H ° = £ S E fimfim + J2 £k C km C km, 
i,m k,m 

n * = E Vk (^ Ckm eikRi + kc -) • 

Here, i is the 4/ site index, k is the conduction electron 
wave vector, and Vu is the hybridization matrix element 
between conduction and localized electrons. £/ and £k, 
both measured from the chemical potential fj,, are the 
excitation energies for 4/ and conduction electrons, re- 
spectively. As a simplification, both types of electrons are 
assumed to have the same angular momentum index m 
with Uf values, m = l...Uf. This will make it possible to 
classify terms in an 1/^/ expansion as is well known from 
the one- impurity model [11]. Finally, the local Coulomb 
repulsion Uf at the 4/ sites has been assumed to be in- 
finitely large so that localized sites can either be empty 
or singly occupied, i.e., the fj m is defined by 



■> i n 



ft II C 1 



J- 



(42) 



A. Perturbation theory 

For the model (41) we first evaluate the effective 
Hamiltonian TL\ in perturbation theory. In contrast to 
the general discussion in Sec. II now we are only inter- 
ested in contributions which renormalize the unperturbed 
part TLq of the Hamiltonian. Using Eq. (13) one finds 



ST~Lo,\ = TLq \ — TLq 



1 E 



(43) 



N 



\s k -s f \>\ 



l^kl 2 

e k - E t 



fim c km, c kmfii 



Note that contributions which involve different / sites 
have been neglected for simplicity. In the following we 
restrict ourselves to the dominant order in Vf so that Eq. 
(43) reads 



8H , 



1 

N 



E 



i,k.m 

£k-£/l>>. 



| ifimfim) 



C lkm C kii 



(44) 



^( C Lri C km)/im// m + ( c ] 1 m Cl[m ) {fimfim) | ■ 

The restriction to contributions of dominant order in vj 
leads to a partial suppression of charge fluctuations on 
/ sites. The contributions, taken into account in Eq. 
(44), describe first the creation of an / and the anni- 
hilation of an c-electron followed by the inverse process. 
Therefore, the rcnormalization contributions in dominant 
order in v$ are based on the existence of a finite prob- 
ability for empty / sites. Furthermore, in Eq. (44) a 
factorization with respect to the full Hamiltonian TL was 
carried out. Note that a factorization with respect to 
TL instead of with TL\ is suggested by the following rea- 
son: As discussed in Sec. II, TL\ results from an unitary 
transformation which is applied to the full Hamiltonian 
TL. Expectation values are not changed by an unitary 
transformation if both operators and the Hamiltonian are 
transformed. Therefore, the factorization with respect to 
the original Hamiltonian TL corresponds to taking into 
account higher-order renormalization. 

Equation (44) can be simplified if we introduce the 
projection operator on an empty / state at site i 



Vo{i) = ~ nL) = fir, 

th 

= 1 — ^2 fimfim 



t 

ill) 



(45) 



(independent of m). Note that the last equation exploits 
the completeness relation on an / site. Here, the / oc- 
cupation number operator, n{ — ^2 m f\ m fi m -, is at the 
same time the projection on the singly occupied / site 
i. According to Eq. (44) the renormalization of the un- 
perturbed Hamiltonian TLo can be written up to second 
order as 



SHo,. 



fc M E fimfim + E fck ' A C km c km + SE X , 



(46) 



Se 



/•A 



fck,A = 



eflet-e/l-A) 



L E 



5E X 



N 



N 



E 



£k-£/ 

l^kl 2 

Ek-Ef 



(Po(i))e(| £k -e/|-A), 
KJ (l-m>(i)» 



xe(| £k - e/ |-A) 



where ds^x and ^£/,a lead to the renormalization of the 
one-particle energies £k.A and £/.a. SE\ is an additional 
energy shift. 
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B. Renormalization equations 

The above result in perturbation can be used to derive 
flow equations for the periodic Anderson model along the 
line discussed in Sec. III. Let us start by formally writing 
down the effective Hamiltonian Ti\ after all excitations 
with energy differences larger than A have been elimi- 
nated. H\ should have the following form 

Wa = Wo,a+Wm, (47) 

T-Co.X = £/,A E fimfim + E £k > X C km C km + #A, 



Wi,A — Pa^i 

= 4= E ^ k (/Lc km e ikR - 



h.c, 



i,k,m 

£ k,A — e/,A I <A 



In the following we are only interested in the renormal- 
ization of the diagonal part of the Hamiltonian. Con- 
sequently, all additional terms which would appear due 
to Eqs. (16) and (17) will be neglected in the following 
for simplicity. Due to renormalization processes the one- 
particle energies £/,a, £k,A, as well as the energy shift E\ 
now depend on the energy cutoff A. The initial conditions 
are those from the original model 

£/,(A=A) = £/, £k,(A=A) = £k, E (\=A) = 0. (48) 

In the following step let us evaluate a new effective Hamil- 
tonian which is obtained after a further elimination of 
excitations within a small energy shell between A — AA 
and A has been done. The new Hamiltonian W(a-aa) is 
obtained from Eq. (46) and possesses the following renor- 
malized parameters 



;/,(a-aa) 



\V k 



N ^ £ k ,A - £/,A 



'km/ J 



(49) 



£k,(A-AA) — £k,A 



1 E' 



Wv 



i?(A-AA)=i?A + ^E' 



£k,A — £/,A 
2 



N 1£ £ k,A-£/,A 

x(l-(V (i)))- 

Here, the prime / above the sums indicates that the re- 
striction (A — AA) < |ek,A — £/,a| < A has to be fulfilled. 

Finally, we would like to perform the stepwise trans- 
formation continuously. By taking the limit AA — > 0, 
we obtain the following flow equations for the parameter 
£/,A, £k,A, and E\ 



de /,A _ »f sr 

k 



E ~ \' K \.. ("km) S ( A - l £ k,A - £/,A| 



dX N^e KX ~e f 



(50) 



de KX _ \V k \ 2 

d\ £k,A — £/J 



(p (i))5(A-| £k ,A-e/,A|), (51) 



dE x 
dX 



= ~ v f 



E 



\V{k)Y 

•k,A — £/,A 



KJ (1 - (V (i))) (52) 



x8(A - |e k ,A - £/,a|)) 
= -N((l-V (i)))^- 

Here, due to the k summation in Eq. (50) it was assumed 
that the renormalization of the one-particle energies Ef t \ 
and £k,A is continuous in contrast to the case of Sec. IV. 
The first two equations are coupled nonlinear differential 
equations for the parameters £f t \ and £k,A- Here, the 
((-function guarantees that only excitation on the energy 
shell contribute. 

To solve Eqs. (50) and (51) first note that due to 
translational invariance the expectation value (Vo(i)) is 
independent of site index i. Therefore, a simple relation 
between Ef.\ and £k,A can be derived by inserting (51) 
into (50) 



CPoii)) 



de 



( /A 



N ^ 

k 



de 



k.A 



dX 



(53) 



Its integration between the lower cutoff A — > and the 
cutoff A of the original model leads to 



(Po(i)) (£/ - if) = E ( n Z™) (ek " £k 



(54) 



where we have defined £/ = E/^a-^o) and ik — £k,(A^o)- 
Note that as before the expectation values in Eq. (54) are 
assumed to be independent of A since they were defined 
with the full Hamiltonian. Eq.(54) is the first equation 
which relates e / and £& ■ To find a second equation let us 
integrate Eq. (51) between the lower cutoff (A — > 0) and 
the cutoff A of the original model. One finds 



£k = £fe 



dX 



'(A^0) £k,A'-£/,A' 

x5(X' - |e k ,A' - £/,A'|) 



(55) 



where Eq. (48) was used. Since the flow of the one- 
particle energies £k,A and £/ ; a was assumed to be conti- 
nous, the main contribution to the integral in Eq. (55) 
should originate from small values of the denominators 
£k.A' ~ £/.A', i-C, from small values of A'. Therefore, 
a reasonable approximation is to replace the quantity 
£k.A' — £/.A' in Eq. (55) by its value at the lower cut- 
off A' « A — > 0. Thus we obtain 



£k 



£k - if 



(56) 
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Note that additional contributions of higher order in 
| Vic | 2 have been neglected. 

Equations (56) and (54) fix the one-particle energies. 
The Hamiltonian 7i(^^ ) reads 

W(A^O) = £/ fLfim + ik C km C km + E (57) 
i.m k,m 



where 



i// ^ |Vfc| (n^ 



Ek = £k 



I^I 2 (W) ) 

£ k - £/ 



£ = TV (l-(P (i)» (£/-£/)• 



(58) 

(59) 
(60) 



As before, the final model Wn_ >q) consists of nonin- 
teracting / and conduction electrons with renormalized 
one-particle energies if and i^. The expectation values 



(Po(i)) 



1 



and (n k ) in Eqs. (58) and (60) can 



again be evaluated from the free energy (for details see 
Appendix B). 

Let us solve Eq. (59) for i^. We find two solutions 



£ k(l,2) 



\ {^k + e f T V / (£k-e7) 2 +4|V/ k | 2 (PoW)| 

(61) 



= 2 i £ k + if T W k } 



which agree with the well-known two heavy quasiparti- 
cle bands known from the slave-boson formalism [12]. In 
the following let us assume that the renormalized / level 
if lies slightly above the chemical potential fi whereas 
the lower quasiparticle band £k(i) is intersected by \i. As 
is easily seen our result is somewhat different from that 
obtained in slave-boson mean-field theory. The renormal- 
ized one-particle excitation e"k in Eq. (61) agrees with the 
lower heavy quasiparticle band of the slave-boson theory 
£k = £k(i)- However, the renormalized / electron exci- 
tation if is different from the upper heavy quasiparticle 
band. Instead we have found a dispersion less excita- 
tion which is located slightly above the Fermi level (sec 
appendix B). Note that for temperature T — both ap- 
proaches lead to same result for thermodynamic quanti- 
ties since at zero temperature only the lower quasiparticle 
part is filled up to the Fermi level. However, static prop- 
erties for T ^ as well as dynamic correlation function 
turn out to be different within the two approaches. 



C. Transformed operators and electron Green's 
function 



To evaluate dynamic and different static quantities we 
have again to apply the transformation on operator quan- 



tities inside the expectation values. For the one-particle 
Green's functions we need 



e Xx cL. e" 



(62) 



and also // m (A). In first order in the hybridization Vk 
one finds for the conduction electron creation operator 



-"km 



T^(QaWi),4 



(63) 



J k?^ 



where ft m is the Fourier transform of the f! m . This re- 

•f Km " ITTl 

suit suggests to make the following ansatz 



C L( A ) = U k,A C km + f k ,A /L 



(64) 



which should be suited to obtain a good approximation 
for cj cm (A). Because cj cm (A), Ck m (A) have to be fulfill the 
Fermi anticummutator relations, one concludes that 

1 = KaI 2 + Ka| 2 (V (i)) (65) 
holds for all values of k. Thereby, the approximation 



ft i 

•> km ! J K 



(n(0> 



(66) 



was used. By using Eq. (64), the conduction electron 
occupation probability reads 



km C km7 



| u k,(A^0)| 



(4m(A-0) C km(A^0)) 



(67) 

(A-0) 



+ K,(A->0)| 



1 



vf + eP s t 



where (...)(^_>o) again denotes the expectation value 
formed with the renormalized Hamiltonian Ti^x~*o) ■ Note 
that 



(fimfim) 



1 



(A-0) 



(68) 



was used in Eq. (67) which differs from the Fermi dis- 
tribution due to definition (42) for the modified creation 
and annihilation operators fj m and /j m . 

On the other hand, the expectation value (c km c k m) 
can also be calculated by functional derivative from the 
free energy (for details see Appendix B). For zero tem- 
perature one obtains 



1 ( g k - if 

2 V VF k 



S(k F - k) 



(69) 



where Wk was defined by Eq. (61). By comparing Eq. 
(67) with Eq. (69) one finds (k<k F ,p^> oo) 
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I |2 1/, £fc - £/ 



and 



|*>k,(A— 0)| 



Gkm(w) 



I 1 + *ZZ1L 



/km 7 



1 ft 



L - (w + irf) 



(70) 




where condition (65) was used. For the transformed / 
electron creation operator one finds 

fL W = (Po(i)> «k,A4 m + u K JL (72) 

where once again Eq. (66) was used. We are now able to 
determine the / electron Green's function 



1 



L (A^o) - (w + it]) 



/L(A-O) 



(73) 



+ / (A^O) 



where new Liouville operators L and with respect 
to H. and Tix were introduced. For the imaginary part of 
the / electron Green's function we obtain 

ImGkm(w) = (Po(i)) {(■Po(i)) \ vk,(\->o)\ 2 $ ~ e k ) 

+ \u k!(x ^ 0) \ 2 S((j-e f )} (74) 

where again k < Uf has to be fulfilled. As is seen from 
Eqs. (70) and (71) for small k values the peak at if 
dominates the imaginary part of the /-electron Green's 
function. The pole at £k becomes important for k values 
near the Fermi momentum kp. 

D. Discussion of the results 

Let us compare our results for the periodic Anderson 
model with those obtained from the slave-boson mean 
field theory. For simplicity we will restrict ourselves to 
the T = case. As already mentioned above, the renor- 
malized energy corresponds to the lower quasiparticlc 
bands from the slave-boson mean field treatment. Fur- 
thermore, the renormalized Hamiltonian TL 



(A-0) 



leads to 



same / state occupation as the slave-boson theory (for 
details see Appendix B). Moreover, the typical Kondo 
peak at the Fermi level is found in the imaginary part 
of the / electron Green's function (74). Therefore, 
the renormalization approach for the periodic Ander- 
son model leads at T = to the same results as the 
slave-boson mean field theory On the other hand, some 
differences occur. The most important one is the addi- 
tional quasi-particlc band at the renormalized / level £/ 
which displaces the upper heavy quasiparticlc band from 
the slave-boson mean-field treatment. Note however that 
in the present approach we have neglected all processes 
which include different / sites. A detailed discussion will 
be given in a forthcoming paper [13]. 

VI. CONCLUSION 



ans. First, the elimination of high energy transition op- 
erators larger than an energy cutoff A leads to a trans- 
formed Hamiltonian 7i,\ which is band diagonal with re- 
spect to the eigenbasis of its unperturbed part. The ap- 
proach, which is based on perturbation theory, can be 
extended to establish a renormalization approach by con- 
tinuously shifting the energy cutoff A to lower and lower 
values. In this way flow equations for the Hamiltonian 
are obtained. The present approach has some similari- 
ties with recently introduced renormalization method by 
Wegner [1] and by Glatzek and Wilson [3,4]. However, 
there are also some substantial differences: The present 
method starts from unitary transformations which are di- 
rectly applied to the Hamiltonian. In contrast, the previ- 
ous approach is formulated in matrix notation. It starts 
from flow equations in differential form for the Hamilto- 
nian as well as for the generator of the unitary transfor- 
mation. The present flow equations for the Hamiltonian 
can be formulated independent from a concrete represen- 
tation. For demonstration the method was first applied 
to the exactly solvable Fano- Anderson model and to the 
Anderson-lattice model. For the first model the exact 
results were found. For the second model some approx- 
imations had to be done. For this case the well-known 
quasiparticle behavior of heavy fermions was rederived. 
In particular, it was shown that there is an almost com- 
plete agreement with the results from the slaved-boson 
mean- field theory [12]. 
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APPENDIX A: FLOW EQUATION IN MATRIX 
NOTATION 



In summary, in this paper we have presented a 
renormalization approach to many-particle Hamiltoni- 



Let us start from the eigenvalue problem of the unper- 
turbed Hamiltonian TLq \ after all transitions with encr- 
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gies larger than A have been integrated out 

H ,a |n(A)) = E& x |n(A)) . (Al) 



The eigenvectors of Ho, a may depend on A. We con- 
sider two degenerate eigenstates \n{\)), \n{\)) belonging 
to the same energy E® A = A . Let us assume that the 
eigenstates have been orthogonalized, i.e., ^n(A)|n(A)) = 
5 n ,h- We multiply Eq. (16) from the left with (n(A)| 
and from the right with |n(A)). For the renormalization 
of Ho, a due to the elimination of excitations within the 
energy shell between A — AA and A one finds 

(n(X)\(H ,(x-AX) -n ,x)\n(X)) = (A2) 
(»(A)|«i,A|m(A)>(m(A)|tt M |n(A)) 



E 

m 

\-AX<\E° X - B m a|^ a 



E° 



E° , 

n, A 



Note that the energy constraint A-AA < \E° X -E^ x \ < 
A may also be written as difference of two 6 functions 



APPENDIX B: CALCULATION OF THE 
OCCUPATION PROBABILITIES 

In this appendix we derive formal expressions for the 
particle number expectation values for the periodic An- 
derson model. Thereby, we restrict ourselves to zero tem- 
perature. Let us assume that the renormalized / level £/ 
lies slightly above the chemical potential /i whereas the 
lower quasiparticlc band £k(i) is intersected by fi. Thus 
for T = 0, according to Eqs. (57) and (60) the ground- 
state energy E g = F{T — > 0) is given by 

E g = E + v f J2 e ~k (Bl) 

= N(1- CP (i))) (e f -e f ) + Vj £k 

k^kp 

where kp denotes the Fermi momentum and £k(i) = £k- 
The ground state energy can be used to calculate the 
zero-temperature occupation probabilities. The func- 
tional derivative of Eq. (57) leads to 



E ••• = (A3) 

\ — A\<\E° A --E", A |<A 

= £••• {6[|< A -< A |-(A-AA)] 
A)} 



-©(Ka-<aI 



\_dEg_ 

N de s 
1 dE a 



1 dE i/f 
Nde] + TV 



"I 



1 dE 
Vf dew 



Eg <** (B2) 

k<k F 



E 

k'<k f 



<9e k 



(B3) 



so that on the rhs of Eq. (A2) the limit AA — > can be 
easily performed. For the lhs it is easy to verify that the 
limit AA — > leads to the derivative of E® A with respect 
to A. Here, the orthogonality of the vectors |n(A)), |n(A)) 
has to be used. Thus, one finds 



dX 



E 



(n(A)|« 1|A |m(A)>(m(A)|« 1)A |n(A)) 



E° - E° 

a m,X 7i, A 



dQ(\E%-El x \-X) 
dX 



In Eq. (A4) the A derivative of the 8 function acts on the 
explicit term linear in A and on the A-dependent energies 
E® x and E^ x . It is obvious that the latter dependence 
leads to higher order contributions in Hi, a- Therefore, 
by restricting one-self to the lowest non vanishing order 
in Hi a one finds 



dE" 



dX 



= -5^(n(A)|Wi lA |m(A)>(m(A)|«i,A|n(A)) 

m 

i 



E° y~E\ 
m,A n , A 



S(\E» X -E»\-X). (A5) 



Note that flow equations for the A dependent matrix ele- 
ments of the interaction Hi, a can be derived in a similar 
way. 



where due to Eq. (60) 



dE 



Bit 



8(r 



Furthermore, by using Eq. (61) one obtains 

^ = 1 (i + £ jL^ll) d lL + l^l 2 d(n{) 
ds } 2 V VF k / de f VF k de f 



(B4) 



(B5) 



(A4) By inserting Eqs. (B4) and (B5) into Eq. (B2) one finds 



= 



dit 



d(n{) 



- VJ - Y fl 

2 N ^ \ 



A-<fe f 



£k ~ Sf 



(B6) 



de f 

We can conclude 

In') = 

[ 1 ' 2 N 



k<kF 



and 



£k ~ £/ 



Vkl 



£/ £f ~ N A* w k 



(B7) 



(B8) 



/->■ <■. 
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Thus the / state occupation of the slave boson mean- 
field theory is reobtaincd [12]. (n^ m ) can be found in the 
same way 



\ n km) 



1 



(B9) 



dij_ 



dew 



N 
v f 

N 



\ £ f 1 



w v 

2 



k'<k f 



By using Eqs. (B7) and (B8) one finally obtains 



Am) =U 1 ~ £ \^- ) 



(BIO) 
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